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Abstract 

We consider a nonlinear convex stochastic homogenization prob- 
lem, in a stationary setting. In practice, the deterministic homoge- 
nized energy density can only be approximated by a random apparent 
energy density, obtained by solving the corrector problem on a trun- 
cated domain. 

We show that the technique of antithetic variables can be used to 
reduce the variance of the computed quantities, and thereby decrease 
the computational cost at equal accuracy. This leads to an efficient 
approach for approximating expectations of the apparent homogenized 
energy density and of related quantities. 

The efficiency of the approach is numerically illustrated on several 
test cases. Some elements of analysis are also provided. 

Keywords: stochastic homogenization, nonhnear problem, variance reduc- 
tion, antithetic variables. 



1 Introduction 



In this article, we consider some theoretical and numerical questions related 
to variance reduction techniques for some nonlinear convex stochastic ho- 
mogenization problems. In short, we show here that a technique based on 
antithetic variables can be used in that context, provide some elements of 
analysis, and demonstrate numerically the efficiency of that approach on 
several test cases. This work is a follow-up of the articles O El [10] where 
the same questions are considered for a linear elliptic equation in divergence 
form. 

The stochastic homogenization problem we consider here writes as follows. 
Let V be an open bounded domain of M'^ and 2 < p < oo. We consider the 
highly oscillatory problem 

mi W (^^,u,Vu{x)^ dx - J f{x)u{x)dx, ueW^'^iV)^ (1) 

for some / and some random smooth field W, which is stationary in a sense 
made precise below, and satisfies some convexity and growth conditions such 
that, for any e > 0, problem ([1]) is well-posed. See Section [LT] below for a 
precise description of the mathematical setting, which has been introduced 
in [m [12]. A classical example that motivated this framework is when 

W{y,u,0 = ^ci{y,u)\^f, 

where a is stationary (see e.g. [HI page 382]). 

In ([1]), e denotes a supposedly small, positive constant that models the 
smallest possible scale present in the problem. For e small, it is extremely 
expensive, in practice, to directly attack ([T]) with a numerical discretization. 
A useful practical approach is to first approximate ([T]) by its associated ho- 
mogenized problem, which reads 

inf (Vm(x)) dx- J f{x)u{x)dx, u G W^'^iV) | , (2) 

and next numerically solve the latter problem. The two-fold advantage of ([2]) 
as compared to ([T]) is that it is deterministic and it does not involve the small 
scale e. 
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This simplification comes at a price. The homogenized energy density W* 
in 02]) is given by an integral involving a so-called corrector function, solution 
to a nonlinear problem (see ([7]) below for a precise formula). As most often in 
stochastic homogenization, this corrector problem is set on the entire space 
W^. In practice, approximations are therefore in order. A standard approach 
(see e.g. [7] in the linear setting) is to generate realizations of the energy 
density W over a finite, supposedly large volume at the microscale, that we 
denote Qn, and approach the homogenized energy density by some empir- 
ical means using approximate correctors computed on Q^. Although the 
exact homogenized density W* is deterministic, its practical approximation 
is random, due to the truncation procedure. It is then natural to generate 
several realizations. However, efficiently averaging over these realizations re- 
quire to understand how variance affects the result. This is the purpose of 
the present article to investigate some questions in this direction, both from 
the theoretical and numerical standpoints. 

Before proceeding and for the sake of consistency, we now present the 
framework of nonlinear stochastic homogenization we adopt, and situate the 
questions under consideration in a more general context. 

1.1 Homogenization theoretical setting 

To begin with, we introduce the basic setting of stochastic homogenization 
we will employ. We refer to [13] for a general, numerically oriented presen- 
tation, and to in m [T7] for classical textbooks. We also refer to [19] and 
the review article [2] (and the extensive bibliography contained therein) for 
a presentation of our particular setting. Throughout this article, {Q, T , P) 



is a probability space and we denote by E(X) = / X{uj)(S^{iJ) the expec- 



tation value of any random variable X G L^(f2,(iP). We next fix (i G N* 
(the ambient physical dimension), and assume that the group (Z'^, +) acts 
on f2. We denote by {jk)k&,<i this action, and assume that it preserves the 
measure P, that is, for all A; G Z'^ and all A^7, P(rfcv4) = P(A). We assume 
that the action r is ergodic, that is, if A G J-" is such that TkA = A for any 
k G Z'^, then F{A) = or 1. In addition, we define the following notion of 
stationarity (see p, Section 2.2]): a function F G Ll^^. (M.^,L^{Q)) is said to 
be stationary if, for all A; G Z"', 

F{y + /c, w) = F{y, Tku) almost everywhere and almost surely. (3) 




3 



In this setting, the ergodic theorem [TSl 1221 El] can be stated as follows: Let 

F G L°° (M."^ , (Q)^ be a stationary random variable in the above sense. For 
k = {ki, k2, . . . , kd) E I^'^ , we set \k\c,o = sup |/ci|. Then 



l<i<d 



T V F{y,TkUj) — >E{F{y,-)) in L^{W^), almost surely. 

+ 1 ) N^oo 



This implies that ( denoting by Q the unit cube in W^) 



F{-,io 



E F{y, ■)dy^ in L°°(M.'^), almost surely. 



The purpose of the above setting is simply to formalize that, even though 
realizations may vary, the function F at point ?/ G M"^ and the function F at 
point y + k, k & Z*^, share the same law. In the homogenization context we 
now turn to, this means that the local, microscopic environment (encoded 
in the energy density W) is everywhere the same on average. From this, 
homogenized, macroscopic properties will follow. 

We now describe more precisely the multiscale random problem ([1]) . The 
domain P is a regular (in the sense its boundaries are Lipschitz-continuous) 
bounded domain of M*^. The right-hand side function / belongs to (P), 
with 1/p + 1/p' = 1 (hence / is indeed in the dual space of L^iV)). For 
any ^ G M'', the random field y,co ^-^ W{y, 0J,C,) is assumed stationary in the 
sense ([3]). We assume that it is continuous (and even C^) with respect to 
the ^ variable, and that it is measurable with respect to the y argument. We 
also assume that there exists C2 > Ci > such that 

Vi/GM^ \fuen, V^GM", c,\^\^<W{y,u,0<C2{l + m- (4) 

Furthermore, we assume henceforth that W is strictly convex with respect to 
the argument ^, in the sense that 

Vr^GM^ V^gM^ ri'^d^W{y,uj,^)r] > a.e. and a.s., (5) 

where d'^W G R'^^'^ is the Hessian matrix of ^ W{y,u,C,). A more de- 
manding assumption is that W is a- convex with respect to the argument ^, 
in the sense that there exists a > such that 

Vt^gM'^, V^gM'^, ri'^d1W{y,u,0r]> a r]^r] a.e. and a.s. (6) 
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Unless otherwise stated, we only assume in the sequel. When needed, we 
will explicitly assume (E]). 

Under (jl]) and (|5]), the variational problem ([1]) is well-posed. In ad- 
dition, the homogenized limit of ([1]) has been identified in [TT| 112] (see 
also [m Theorem 3.1]): the unique solution u'^{-,ijj) to ([1]) converges (weakly 
in W^'^{V) and strongly in Lp(T'), almost surely) to some deterministic func- 
tion u* G ly^'^(P), solution to ([2]), where the homogenized energy density 
W* is given, for any ^ G M'^, by 

W*iO = lim inf / W iy,co,^ + Vwiy)) dy, w G Wj^^iqA 

N^cc [\Qn\Jqj, ^ J 

(7) 

where Qn = {—N,Ny and where W^^{Qn) denotes the set of functions 
that belong to W^^^(R'^) and are QAf-periodic. The convergence in ([7]) holds 
almost surely. 

1.2 The questions we consider 

In practice, we cannot compute W*{C,), and have to restrict ourselves to finite 
size domains. We therefore introduce 

W^^(a;,0:=mf(^ / W iy,oj,^ + Vw{y)) dy, «; G l^i'^(g^)) (8) 

and readily see from (I7j) that 

W^'^(0= lim W;,iu,0 a.s. 

As briefiy explained above, although W* itself is a deterministic object, its 
practical approximation is random. It is only in the limit of infinitely 
large domains Q n that the deterministic value is attained. This is a standard 
situation in stochastic homogenization. 

Many studies have been recently devoted (at least in the linear case) 
to establishing sharp estimates on the convergence of the random apparent 
homogenized quantities (computed on Qn) to the exact deterministic ho- 
mogenized quantities. We refer e.g. to [TJ [15] and to the comprehensive 
discussion of [HI Section 1.2]. We take here the problem from a slightly 
different perspective. We observe that the error 
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is the sum of a systematic error (the first term in the above right-hand side) 
and of a statistical error (the second term in the above right-hand side). 
We focus here on the statistical error, and propose approaches to reduce the 
confidence interval of empirical means approximating E [W^jv(-, 0] (o^ similar 
quantities), for a given truncated domain Qn- 

Recall that a standard technique to compute an approximation of E [W^(-, ^)] 
is to consider several independent and identically distributed realizations 
of the energy density W, solve for each of them the corrector problem ([8]) 
(thereby obtaining several i.i.d. values Wp"^{u,C,)), and proceed following a 
Monte Carlo approach: 

^ 2M 
m=l 

In view of the Central Limit Theorem, we know that our quantity of interest 
E [M/Ar(-,0] li^s in the confidence interval 

hu - 1-96 -== , I2M + 1-96 -== 

with a probability equal to 95 %. 

In this article, we show that, using a well known variance reduction tech- 
nique, the technique of antithetic variables [2T1 page 27], we can design a 
practical approach that, for finite and any vector ^, allows to compute a 
better approximation of E[P1/'^(-,^)] (and likewise for similar homogenized 
quantities). Otherwise stated, for an equal computational cost, the approach 
provides a more accurate (i.e. with a smaller confidence interval) approx- 
imation. We thereby extend to this nonlinear convex setting the results 
of El [ID] obtained in the linear case. 

Our article is articulated as follows. In Section 12. we describe the pro- 
posed approach, and state our main results. The ingredients to prove these 
results are collected in Sections 12. 2[ 12.31 and 12.41 The actual proof of our 
main results is performed in Section 12.51 We make there several structural 
assumptions on the form of the energy density W to obtain these variance 
reduction results. In Section 12. 6[ we describe a general class of examples for 
which our assumptions are indeed satisfied. We next turn in Section [3] to 
some illustrative numerical examples, where we demonstrate the efficiency of 
the approach, even in cases where the theoretical analysis is incomplete. 
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2 Description of the proposed approach and 
main results 



2.1 Statement of our main results 

This section is devoted to the presentation and the analysis of our approach. 
We first focus on estimating the expectation E ^)] of the apparent ho- 

mogenized energy density (see Section [2XT]). Our variance reduction result, 
Proposition [H shows that the technique of antithetic variables is indeed ef- 
ficient. As often the case, it is difficult to quantitatively assess how efficient 
the approach is, and this will be the purpose of the numerical tests described 
in Section [3] to address this question. 

We then turn to the estimation of the first (and next second) derivatives of 
W^{-,^) with respect to ^. These quantities naturally appear when one solves 
the convex homogenized problem ([2]) (approximating W* by ^^^(0;, ■)), e.g. 
using a Newton algorithm. For these two quantities, our result is restricted 
to the one-dimensional setting. See Section 12.1.21 and Proposition [2] for the 
first derivative, and Section [2.1.3l and Proposition [3] for the second derivative. 

Sections 12.2112. 31 12^ and l2.5l are devoted to the proof of the results stated 
here. In Section 12.61 we discuss an explicit class of energy densities W that 
falls into our framework. 

2.1.1 Variance reduction on the homogenized energy density 

In this section, we make the following two structure assumptions on the 
rapidly oscillating field of ([T]). First, we assume that, for any A^, there 
exists an integer n (possibly n = \Qn\, but not necessarily) and a function 
A, defined on Qn x R" x R'', such that the field W{y,u,C,) writes 

WyeQN, V^gR^ W{y,u,0 = A{y,X,{u),...,Xn{co),0 a.s., (9) 

where {A^A;(<^)}i<a;<„ are independent scalar random variables, which are all 
distributed according to the uniform law W[0, 1]. In general, the function 
A, as well as the number n of independent, identically distributed variables 
involved in depend on A^, the size of Qn, although this dependency is 
not made explicit in 

Second, we assume that the function ^ in ([9]) is such that, for all ?/ G Qat 
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and all ^ G M'^, the map 



(10) 

is non- decreasing with respect to each of its arguments. 

Proposition 1. We assume flUl) - ffTU]) . Let W^{u,C,) be the approximated 
homogenized energy density field defined by ([HD- We define on Qn the field 

W^^\y, u, := A{y, 1 - X,{u), . . . , 1 - X„(a;), 0, 

antithetic to W defined by (Q. We associate to this field the approximate 
homogenized energy density field W^^''^{u,^), defined by ([HD (replacing W by 
Set 

wu^, ■■= \ {w'n{^. + wT^\uj, 0) ■ (11) 

Then, for any ^ G M'^, 



E 



=E[iy;^(-,0] and W;,{-,0 <^Var[iy*(-,0]. (12) 



Otherwise stated, W^iuj,^) is a random variable which has the same expec- 
tation as W^{(jj,^), and its variance is smaller than half of that ofW^{u,^). 

As mentioned above, this result generalizes [6i Proposition 2.1] to the 
nonlinear convex variational setting considered here. 

Before proceeding, we briefly explain the usefulness of the above result 
for variance reduction techniques. Assume we want to compute the expec- 
tation of W^{(xJ,^), for some fixed vector ^ G M'^. Following the classical 
Monte-Carlo method recalled in Section [L2| we estimate E[iy^(-,^)] by its 
empirical mean. To this end, we consider 2M independent, identically dis- 
tributed copies {Wm(|/, 1^, 0}i<m<2M raudom field W{y,u,^) on Q^. 
To each copy Wm, we associate an approximate homogenized energy density 
W^'^{u,^), defined by ([H]). We next introduce the empirical mean 



2M 

i 

'2M 



^E^r(c^,0, (13) 

m=l 

and consider that, in practice, the mean E ^)] is equal to the estimator 

I2M within an approximate margin of error 1.96-^^ ; . 

V2M 
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Alternate to considering (fT^ . we may consider 

1 _ 

hM = j^j2^n^,o, (14) 

m=l 

where W^"^ is defined by (fTTl) . Again, in practice, the mean E [W^Ar(-,0] = 



E 



is equal to I2M within an approximate margin of error 1.96 



Var 



Observe now that both estimators f|T3|) and f|T4l) are of equal cost, since they 
require the same number 2M of corrector problems to be solved. The accu- 
racy of the latter is better if and only if Var W^{-,^) < -Var [W^{-,^)], 
which is exactly the bound (|T2|) of Proposition [TJ 



2.1.2 Variance reduction on the first derivative of the homoge- 
nized energy density 

Restricting ourselves to the one-dimensional setting, we now state a variance 
reduction result for the estimation of K[^d^W^{-,^)]. Note that, to distin- 
guish derivatives with respect to y from derivatives with respect to C,, we keep 
the notation d^W, even though we are in the one-dimensional situation. 

We again make the structure assumption ([9]), and observe that it implies 
that 

\/yei~N,N), V^gM, ^d^W{y,u,0 = Ai{y,X,{co),...,X^{u),0 a.s., 

where {-^fc(i^)}i<fc<„ are scalar i.i.d. random variables, which are all dis- 
tributed according to the uniform law W[0,1], and where the function Ai, 
defined on (— A^, A^) x M" x M, is given by 

A,{y,x,0 = ^d^A{y,x,0. (15) 

In addition, we assume that, for all y G (— A^, A^) and all G M, the map 

(xi, . . . , x„) G M" H> Ai{y, xi, . . . , Xn, (16) 

is non- decreasing with respect to each of its arguments. 

We recall that the function ^ (-)■ W{y, uj,^) is strictly convex (see assump- 
tion (^) and satisfies (jl]). It therefore has a unique minimizer ^o{y, u). In the 
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sequel, we consider energy densities such that this minimizer is independent 
of y and uj. Without loss of generality, we can assume that = 0. We thus 
consider energy densities W such that 

^ I—)- W{y,uj,E,) attains its minimum at ,^ = 0, a.e. and a.s. (17) 

Proposition 2. Let d = 1, and assume (jH]), (fT^ . f lTB]) and f lT7|) . We intro- 
duce 

:= I i^d^W^'* {u, + ^d^W;, {CO, 0) , (18) 

where W^^'* (uj,^) and {^,0 are defined as in PropositionUi Then, for 
any ^ G M, 



E 



=E[eW(-,0] and Y^T ^d^W;,{;0 < ^Var W(-, 0] 

(19) 



2.1.3 Variance reduction on the second derivative of the homog- 
enized energy density 

Considering again the one-dimensional setting as in Section I2.1.2[ we even- 
tually state a variance reduction result for the estimation of E [c?|W^Ar(-, i^)] ■ 
Recall that, for any y and u, the map ^ (-)■ d^W{y, w, ^) is increasing. We 
can therefore introduce its reciprocal function ( h-?- tlj{y,uj,(), which is also 
increasing. 

We again make the structure assumption and observe that it implies 
that, for any y G (— A^, A^) and any C ^ 

dfW {y,u,^p{y,u,C)) = ■A2iy,Xi{u), . . . ,X„(w),C) a.s., 

where {-^fc(w)}]^<;j<„ are scalar i.i.d. random variables, which are all dis- 
tributed according to the uniform law W[0, 1], and where the function A2, 
defined on (— A^, A^) x M" x M, is given by 

A2{y, X, C) = diA {y, x, [d^A{y, x, ■)]"' (0) , (20) 

where ( ^ [d^A{y,x, {() is the reciprocal function of ^ h-j- d^A{y,x,^). 
In addition, we assume that, for all y G (— A^, A^) and all C ^ map 

{Xi,...,Xn) GR"h^^2(y (21) 

is non- decreasing with respect to each of its arguments. 
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Proposition 3. Let d = 1, and assume ©, (US}, (USD, (EHD one? ([21]). We 
also assume that (1171) holds, and that 



^ I— 7- d'^W {y, u,^) is non decreasing for ^ > 
and non increasing for < 0, a.e. and a.s. 



(22) 



We introduce 



where W^*'* '^'^^ '^n i^^O ^'^^ defined as in PropositionUl Then, for 

any ^ G M, 



E 



dlW*Ni-^0\ = E and Var < 2^ar 

(23) 



The density W{y,u,^) = a(?/,a;)|^|^, where a is positive and bounded 
away from zero and p > 2, typically satisfies the assumption (l22l) . 



2.2 Classical results on antithetic variables 

We first recall the following lemma, and provide its proof for consistency. 
This result is crucial for our proof of variance reduction using the technique 
of antithetic variables, performed in Section 12.51 

Lemma 4 ([21], page 27). Let f and g be two real-valued functions defined 
on M", which are non- decreasing with respect to each of their arguments. 
Consider X = {Xi, . . . ,X„) a vector of random variables, which are all in- 
dependent from one another. Then 

Cov(/(X),(?(X))>0. (24) 

Proof. This lemma is proved by induction. We treat the one- dimensional 
case (n = 1) below, and we refer to [6l Proof of Lemma 2.1] for the induc- 
tion. Consider X and Y two independent scalar random variables, identically 
distributed. Both functions / and g are non-decreasing, so 

ifiX)-f{Y)) {g{X)-g{Y))>0. 
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We now take the expectation of the above inequahty: 



E(/(X) g{X)) + E(/(F) g{Y)) > E(/(F) g{X)) + E(/(X) g{Y)). 



As X and y share the same law, and are independent, this yields 



E(/(X) g{X)) > E(/(X)) Eig{X)) 



and fl24p follows for n = 1. 



□ 



The following result is a simple consequence of the above lemma (see 
e.g. |6] for a proof). 

Corollary 5 ([21]). Let f be a function defined on M", which is non- decreasing 
with respect to each of its arguments. Consider X = {Xi, . . . ,X„) a vector 
of random variables, which are all independent from one another, and dis- 
tributed according to the uniform law U[0, 1]. Then 



where we denote 1 - X = (1 - Xi, . . . , 1 - X„) G MJ". 

Proof. Choosing g{xi, . . . , x„) = —/(I — Xi, . . . , 1 — a;„) in Lemma HI we 
obtain that 



Cov(/(X), /(I - X)) = Cov(/(Xi, . . . , X„), /(I - Xi, . . . , 1 - X„)) < 0. 




We next observe that 



VarQ (/(X) + /(l-X)) 



= ivar(/(X)) + icov(/(X),/(l-X)) 
< ^Var(/(X)), 



where we have used that Var(/(X)) = Var(/(1 — X)). 



□ 
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2.3 Derivatives of the corrector and of the homoge- 
nized energy density 

We now introduce the correctors as the solutions to ([H]): 

w;^(-,a;,0:=arginf(/ W{;u,^ + Vv), veW^^iQ^), [ v = o] . 

In this section, we derive some useful expressions for the derivatives with 
respect to ^ of and of W^. 

The first order optimahty condition in (|H]) reads 

yheW^^iQN), [ {Vhfd^W{-,co,^ + Vw'')=0. (25) 
We deduce from that condition that 

W(a;,0 = ^/ d^W{;u,^ + Vw''), (26) 

and we note that we do not need to know d^w^ to compute d^W^. Comput- 
ing the derivative of this equahty with respect to S,-, we obtain that 

5|W^^(c^,0 = ^ / (ld + 9gVw^)9|iy (■,a;,e + Vw;^) (27) 



with the convention that r^fVw^l ., = — - — for I < j,k < d. We can 

actually obtain a somewhat more symmetric expression. Computing the 
derivative of fl25|) with respect to ^, we indeed see that 

-iheW^/iQN), [ {ld + d^Vw^)dlW {■,uj,^ + Vw^)Vh = 0. (28) 

JQn 

We then infer from (1271) and (EHD that 



dlW]^{u,0 = T^^ f {ld + d^Vw'')dlW{;u,^ + Vw''){ld + d^Vw''f 
\^n\ Jqj^ 

(29) 
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Remark 6. Using the same kind of arguments, we see that the function 



gj = G W^^{Qn) is solution to the variational formulation 
\/h G W^^'iQN), [ {Vhf dlW (-, io,^ + Vw"") Vg, 

JQn 

Suppose that W is a-convex (i.e. satisfies IQ). Then problem (!30l) is well- 
posed and allows to uniquely determine (up to an additive constant) gj, by 
solving a linear elliptic partial differential equation. 

Combined with fl29l) . this remark provides a practical way to compute 
d'^W^{u,C,) without using any finite difference approximation in ^. 

We finally note that, in view of ( l26l) . we have 

^.d^WU^,0 = j^ [ ^■d^W{;co,^ + Vw'^). (31) 
Likewise, in view of ( 129|1 . we see that 

^ / [e + v(e-9,^^)]^9|iy(-,a;,e + v^^)[e + v(e-9,t.^)]. 

\Qn\ Jq^ 

(32) 



2.4 Monotonicity properties 

Our goal in this section is to establish monotonicity properties for the ho- 
mogenization process. Such properties are indeed useful to apply Corollary [5] 
and therefore prove variance reduction. 

To simplify the notation, we assume in this section that we are in a 
periodic setting. For any ^ G M"^, the function y i— )■ W{y,^) is supposed to 
be Q-periodic (with Q = (0,1)^^), to satisfy the growth condition (jlj) and 
to be strictly convex with respect to ^. The associated homogenized energy 
density is then given by 

W%0 = ^^iy^W{y,^ + Vw{y))dy, weW^^iQ), ^^ = o}- (33) 
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We first show a monotonicity property on the homogenized energy density 
in Section 12.4. 1[ Next, restricting ourselves to the one-dimensional setting, 
we show monotonicity properties for the first and the second derivative of 
the homogenized energy density (see respectively Sections 12.4.21 and I2.4.3p . 

2.4.1 On the homogenized energy density 

The following result is an extension to the nonlinear setting of a well-known 
result in the linear setting (see [231 page 12]). 

Lemma 7. Suppose that the fields Wi and W2 satisfy 

Ve G W2{y,i) > W^i(y,0 a.e. on Q. (34) 

We denote Wi and W2 the corresponding homogenized energy densities, de- 
fined by ( 15^ . We then have 

Ve G W^iO > W^iO- (35) 

Proof. Fix ^ G M.'^. For any v G W^^{Q) with / t> = 0, we have that 

Jq 

W^iO < I ^1 (y, e + ^v{y)) dy< [ W2 {y, e + Vv{y)) dy. 
Jq Jq 

Taking the infimum over v, we obtain the claimed result. □ 

Remark 8. Consider the case of an energy density that is positively homoge- 
neous of degree p with respect to its variable ^, that is such that W{y, A^) = 
IM'^W^y,^) for any y G M.'^, ^ G M*^ and A G M. A typical example is 

= ~(^{y)\^\^- '^^^'^ have, for any y and ^, that 
P 

^■d^W{y,0=pW{y,0 and edlW{y,0^ = p{p - l)W{y,0. (36) 
Using successively (13T1) . ( l25i) and ( |36i) . we obtain that 

Jq 

= [ {^ + Vw)-d^W{-,^ + Vw) 
Jq 

= P [ W{;^ + VW) 

Jq 

= pW^iO, (37) 
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where w is the corrector, solution to f l33p . 

We next observe that, for any \ eM, we have w{-, A^) = \w{-,^). Thus, 
for any y, the map ^ i— )■ w{y,^) is homogeneous of degree one, and therefore 
^ ■ d^w = w. We thus infer from fl32l) . using fl36l) . that 



edlW^m = [[^ + VwfdlW{;^ + Vw)[^ + Vw] 

JQ 

= P(P-I) [ W{-,^ + Vw) 
Jq 

= p{p-l)W\0- (38) 

Consider now two fields Wi and W2 that are positively homogeneous of de- 
gree p with respect to the variable ^ and satisfy fl51|) . Then we deduce 
from dSSD, dSZD dMD that, for all ^ G M"', 

^■d^W^{0>^-d^W^{0 and fdlW^m > fd^W^OC 

2.4.2 On the first derivative of the homogenized energy density 

We now establish a monotonicity result on the derivative of W*{C,), in the 
one-dimensional setting. 

As in Section [2.1.21 (see flTTl) ). we consider energy densities W such that 

C, i-> W{y,C,) attains its minimum at ,^ = for almost all y E Q. (39) 

Lemma 9. Let d = 1, and consider two energy densities Wi and W2 satis- 
fying fl39l) . and such that 

Ve e R, ^d^W2{y,0 > ^d^W^{y,0 a.e. on (0, 1). (40) 

We denote Wi and W2 the corresponding homogenized energy densities, de- 
fined by (155]) . We then have 

veeM, eW(0>e5?W^r(0- (41) 

Proof. We first claim that 

d(W*{^) has the same sign as ^. (42) 
To prove this, we note that the corrector equation reads (see fl25l) ) 
d 

dy 
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d^W(y,^ + ^{y,0 



on (0, 1), ti'(-, ,^) is 1-periodic. 



We therefore see that d^W iy,C, + -i— ) is independent of y, and us- 



Let ^ I—)- ip{y,^) be the reciprocal function of ,^ i— )■ d^W{y,^), which exists 
and is increasing thanks to the strict convexity of ^ ^ W{y,^). We deduce 
from the above equation, after integration over (0, 1), that 



We are now in position to prove ( H2|) . Indeed, we first note that (!39|) . that 
reads di:W{y,^ = 0) = 0, imphes that ^lj{y,0) = 0. If d^W*{0 > 0, then 
ip^y, d^W*{^) > ipi^y, 0) = 0, hence, integrating over (0, 1) and using ( H3ll . we 
obtain ^ > 0. Likewise, di:W*{^) < implies that ^ < 0. The claim is 
proved. 

To proceed, we see that the assumption (HOj) equivalently reads, using the 
reciprocal functions. 



We now prove (14T|) by contradiction. Assume that ^d^W2{^) < ^<9gW^i*(0 
for some ^ G M. Without loss of generality, we can assume that ^ > 0, 
and therefore 5gM^2*(0 < Using (H^ . we additionally have < 

^?^2*(0- Using that ( ^ ^2(2/, C) is increasing and fj4i|) with ( = d^W^{^) > 
0, we have 



ing ( 126|) . we obtain that 





(43) 



vc e R, (My, C) < (My, C) a.e. on (0, 1). 



(44) 



V^2 {y, d^w^iO) < ^2 {y, d^w^iO) < ^1 (y, d^w*{0) ■ 



Integrating over (0, 1) and using (H3l) yields 




i^2{y,d^W^iO)dy< [ i^,{y,d^W^iO)dy = ^, 



and we reach a contradiction. This concludes the proof. 



□ 
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2.4.3 On the second derivative of the homogenized energy density 

We next turn to monotonicity properties of the second derivative of the 
homogenized energy density. As in Section [2.4.2[ we consider energy densities 
satisfying f p9|) . We additionally request that, almost everywhere in (0, 1), 

^ I— )■ d'^W (y, is non decreasing for ^ > , . 

and non increasing for ^ < 0. 

Lemma 10. Let d = 1, and consider two energy densities Wi and W2 satis- 
fying 1^, (HOD, (US} and such that 

VCgM, dlW2{y,MyX))>dlW^{y,MyX)) a.e. on (0,1). (46) 

We denote and W2 the corresponding homogenized energy densities, de- 
fined by 032]) ■ We then have 

G M, dlW'^{i) > dlW^. (47) 
We recall that C ^ "^{u, C) is the reciprocal function of ^ i-)- d^W{y,^). 
Proof. We first compute the derivative of (H3|) and obtain 

1 = f dy 

dlW^iO Jo dlW[y,ij{y,d^W*m' ^ ^ 

It is sufficient to prove (l47j) for ^ > 0. Using (HTl) and the fact that ipi and 
d'^Wi are non-decreasing with respect to their second argument, we have 

dlW, {y, ^1 {y, d^W^m < dlW, {y, ^1 {y, d^W^m ■ 

Using (l46l) for ( = d^W2{^), we deduce that 

dlW, {y, ^1 {y, W(0)) < dlW2 {y, ^2 {y, W(0)) ■ 



In view of f l48|) . this inequality readily implies fl47j) for ^ > 0. This concludes 
the proof □ 

2.5 Proof of Propositions [1], [2] and [3] 

Now that we have collected all the necessary ingredients, we are in position 
to prove our main results. 
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2.5.1 Variance reduction on the homogenized energy density 

Proof of PropositionUi As 1 — Xk{uj) and Xk{uj) share the same law, so do 
the fields W and W^^^ on Qn- Hence, the homogenized fields W^{u,C,) and 
W^^'*{uj,C,) share the same law, and we obtain the first assertion of f|T2l) . 

We now choose a vector ^ e R'^, and denote by Vff the operator that 
associates to a given Q^r-periodic energy density the homogenized energy 
density evaluted at ^. We see from ([8]) that W^{u,^) is the effective energy 
density (evaluated at ^) obtained by periodic homogenization of W\y^Qj^: 

W;,iuj,0 = K[Wi;UJ,\^Q,] a.s. (49) 
Using the function ^ of ([9]), we introduce the map 

/ : R" ^ R 

see that f{X{u)) = W^{u, ^) and that, using the definition ffTTj) of W^{uj, C,), 
we have 

i ifiXico)) + /(I - Xico))) = \ [W^^iuj, + WT^\uj, 0) = W^n{^, i). 

(50) 

We now infer from Assumption ffTOj) that, for any y E Qn and any C G R*^, 
the function A{y, ■, C) is non- decreasing with respect to each of its arguments. 
In view of Lemma [TJ we obtain that / is non- decreasing. 

We are thus in position to use Corollary 0, which yields 
Var Q (/(X) + /(I - X))] < Ivar (/(X)) . 



Using ( 15011 . we obtain 
Var (^W7;,(-,0)=Var 



i(/(X) + /(l-X)) 



<ivar(/(X)) = ivar(W^;y(-,0) 



which concludes the proof of the second assertion of (fT2!) and of the Propo- 
sition [U □ 

Remark 11. Following Remarkl^ consider a positively homogeneous energy 
density W . We have shown there that ■ d^W^^u,^) and ^'^d'^W^{u,^)^ are 
equal (up to a multiplicative constant) to W^{u,^). Thus, under Assump- 
tions fl9|) -f lTOj) ■ variance reduction holds for these two outputs as well. 
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2.5.2 Variance reduction on the first derivative of the homoge- 
nized energy density 

Proof of Proposition The proof follows the same lines as that of Proposi- 
tion [H 

As 1 — Xk{u) and Xk^u) share the same law, so do the fields W and 
VF^°* on Qn. Hence, the quantities ^d^W^jico,^) and ^d^W^""^'* {to , ^) share 
the same law, which implies the first assertion of (fTOl) . 

To prove the second assertion, we again make use, as in the proof of 
Proposition [H of the operator Vf^ that associates to a given QTv-periodic 
energy density the homogenized energy density evaluted at ^ (here, Qn = 
{—N,N)). Expression holds. Choosing a vector ^ G M, we introduce 
the function 

/ : M" R 

which obviously satisfies /(X(w)) = ^d^W^^u,^). Using the definition (fT8l) 
of C^^^^y have 

^ [/(XM) + /(I - x{u))] = \ [id^w;,{uj,i) + id^wT'%u:,i)\ = 

(51) 

We now infer from Assumption ( |T6|) that, for any y G (—X, X) and any 
^ G M, the function Ai{y, ■,() is non-decreasing with respect to each of its 
arguments. In view of Lemma [H we thus obtain that / is non-decreasing. 

Using CoroUaryia we write that Var Q (/(X) + /(I - X))^ < ^Var (/(X)). 
In view of (ISTI) . we recast this inequality as 

Var 



eW(-,0] <^Var(eW(-,0) 



and therefore obtain the second assertion of ffTI?]) . This concludes the proof 
of Proposition [2l □ 



2.5.3 Variance reduction on the second derivative of the homog- 
enized energy density 

Proof of Proposition 0. The proof follows the same lines as the proof of 
Proposition O Using Assumptions ( 1T6|) and ( l2Ti) . we see that Assump- 
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tions fl40l) and fl46p of Lemma [TU] are satisfied. The monotonicity result 
of Lemma [TU] next allows to use Corollary El which implies f l23p . □ 

2.6 Examples satisfying our structure assumptions 

Before proceeding to the numerical tests, we give here some specific examples 
of fields W that satisfy the above assumptions. We consider the case 

W{y,u,0 = a{y,ujy-^ + ciy,uy-^, p>2, (52) 
p 2 

with c{y,u) > and a{y,u) > a_ > a.e. and a.s., and provide suf- 
ficient conditions on the scalar fields a and c for the structure assump- 
tions ([9D, (IIOD, (IISD and (EID to be satisfied. Note that (117]) and ([22]) are 
already fullfilled. 

Consider two families (afc(u;))^g^d and {ck{co)) f.^^d of independent, iden- 
tically distributed random variables, and assume that 

a{y,uj) = lQ+k{y)ak{uj), c{y,uj) = ^ lQ+k{y)ck{uj), (53) 
keZ'i fcez<* 

where Q = (0, l)'^ and Q + k is the cube Q translated by the vector k G Z^. 
The scalar field a{y, u) is therefore constant in each cube Q + k with i.i.d. 
values afc(a;), and likewise for c{y,uj). 

We assume that there exist a > and (3 < oo such that, for all k G Z*^, 
< a < ttkioj) < /3 < +00 and < Ckioj) < /3 < +oo almost surely. 
Consequently, ([!]) holds. 

Introduce now the cumulative distribution functions Pa{x) = z/a(— oo,x), 
where z/^ is the common probability measure of all the a^, and next the non- 
decreasing functions fa{x) = inf {z;Pa{x) > z}. Then, for any random vari- 
able X°'{u) uniformly distributed in [0, 1], the random variable fa{X°-{u)) is 
distributed according to the measure Ua- As a consequence, we can recast ([53]) 
in the form 

a(y,a;)= 5^lQ+fc(y)/„(X»(a;)), 

where {X^{u))j^^^d is a family of independent random variables that are all 
uniformly distributed in [0,1], and /„ is non- decreasing. We can proceed 
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likewise for the variables c^. This yields an example where (Q, fllUI) and fll6l) 
hold. In particular, the function ^ of ([H]) reads 



where In = {k & s.t. Q + A; C Qn} and = {a^fclj^g/^- As shown in [B], 
more general fields a{y,u) (where random variables may be correlated) also 
fall into this framework. 

In what follows, we prove that, under assumptions (!52|) and (153|) . and if 
p < 3, the structure assumption (!2T|) holds. Without loss of generality, we 
may assume that y G (0, 1), and write that 

IflP |£|2 

Vi/G(0,l), ^(a;„,Xe,0 = a— + c^, 

p z 

with a = faix^) and c = /c(a;3)- By a slight abuse of notation, we keep 
implicit the dependency with respect to y, work with d and c rather than Xa 
and Xc, and write 

^(a,c,0 = a— + c^- 
p 2 

We compute 

and denote C f^l*^, c, C) the reciprocal to the function ^ i-t- (9^.4(0, c,^): 

C = d \g{d,c,C)f~^ 9{a,c,C) + cg{d,cX)■ 
Th.e function A2 of (!2T|) then reads 

^2(0, c, C) = (p - l)a Is-la, c, C) + c. 

We are left with showing that A2 is non- decreasing with respect to d and c. 

A first remark is that since g{d,c,() has the same sign as ( (recall that 
a > and c > 0), we may as well restrict ourselves to C > and g{d, c, () > 0. 
We hence have 

^2(0, c,C) = ip-l)dg{d,c,CY~'^ + c, 

C = dg{d,c,(y~^ + cg{d,c,C)- (54) 
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We first compute the derivative of A2 witli respect to a: 

^ = (p - l)g{-a, c, C)^-^ + {p-l){p- 2)a^?(a, c, CY''^- 

Og 

Using (|5^ to compute — , we obtain tliat 

da 

dA.2 

and since p > 1, a > and g > 0, we deduce tliat —;— > 0. 

oa 

We next compute tfie derivative of A2 witli respect to c. Using again flMl) 

dg 1 . , 

to compute — , we obtam tiiat 
oc 

{c+ip- l)-ag^~') ^ = c - (J9 - l)(p - 3)ag^~\ 

Recall that a>0, c>0, p>l and g > 0. We have assumed that p < 3, 

dA.2 

and therefore deduce from the above relation that > 0. The structure 

oc 

assumption ( 12T|) hence holds in that case. 

Remark 12. The argument above also shows that the case 

Wiy,u,0 = aiy,u)^-^, 
P 

along with assumption fl5^ . falls in our framework, for any p > 2. 
It is likely that other settings, such as 

Wiy, u, = (aiy, u) + c{y, u;)) ^ + c{y, a;) 

p I 

along with assumption (!53|) . where at and Ck are all independent random 
variables, also fall in our framework. We will not pursue in this direction 
here. 



3 Numerical results 



Our numerical experiments are presented in Section 13. 2[ and discussed in 
details in the subsequent sections. In Section 13. ![ we first discuss the algo- 
rithm we used to solve the variational problem (|8]) that defines the apparent 
homogenized energy density. 
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3.1 Newton algorithm to solve the truncated corrector 
problem 

As mentioned above, the corrector problem ([HD is a convex minimization 
problem, which has been well studied in the literature (see e.g. [21 El [ISl I2D] ) ■ 
We explain here how we proceed in practice to solve this problem, assuming 
that W is not only strictly convex, but actually a-convex (i.e. satisfies (Q). 

To simplify our exposition, we use the notation of the Q-periodic case, 
where the corrector problem is (1331) . We introduce some basis functions 
{v^ijjg/ (e.g. finite element functions) where ifi G W^^{Q), and the finite 
dimensional space Vh = Span {ipi,i G /}. Consider the functional 

J(w) = J {{w,\^,) = [ Wiy,^ + Vwiy))dy 

JQ 

defined on V/j, with 

i&I 

and the variational problem 

inf J{vh). (55) 

This problem has a unique solution (denoted Wh G Vh) up to the addition of 
a constant. The quantity is well-defined, and is the finite-dimensional 
approximation of Vw, where w is the solution to (!33|) . 

In practice, problem (!55|) is solved using a Newton algorithm. We see 
that 

dJ d'^J 

{w) = D^{ipj) and - — - — [w) = H^{ipj,ipk) 



where 



dwj dwjdwk 

DM= [ Vip{y)-d^W{y,^ + Vw{y))dy 

JQ 



HUv,^)= I {Vv{y)fdlW{y,^ + Vw{y))VtP{y)dy. 



and 

JQ 

The Newton algorithm consists in defining w™'*'^ G Vh from G Vh by the 
following linear elliptic problem: find w™^^ G Vh such that 

v^^ G Vh. - <,^) = -D.^^{e). 
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Again, w"^^ is uniquely defined up to the addition of a constant. 

Tlie finite-dimensional problem f lSSp is a-convex, and W is smooth with 
respect to ^: the Newton algorithm hence locally converges (quadratically), 
and lim Vif™ = Vwh- 

In practice, we consider a sequence Th of meshes on Q, and set Vh = 
Pjj(Q) = {vh € C{Q) s.t. VT G Th,Vh is affine on T}. By classical finite el- 
ement results, we know that lim ||Vwh — Vw||lp(q) = (see e.g. [3] and 

also [251 [1]). 

3.2 Overview of numerical results 

We have considered three test-cases of the form flS21) - flS51) . namely 

IflP |f|2 

W{y, 00, = a{y, u;)^ + c{y, a;)^^ 
p 2 

with a{y,u) = ^ lQ+fc(?/)afc(u;) and c{y,u) = ^ lQ+fc(y)cfc(cu), 

with p = 4, in dimension d = 2. The random variables follow a Bernoulli 
distribution: P(afc = a) = P(afc = /9) = 1/2, with a = 3 and /3 = 23. The 
value of the field c is chosen as follows: 

• Test Case 1: in this first test case, c{y,u) = 0. The problem is thus 
strictly convex but not a-convex. In addition, the energy density is 
positively homogeneous of degree p, hence Remarks [H and [TT] apply. 

• Test Case 2: the second test case corresponds to c{y, u) = 1. The prob- 
lem is then a-convex, and highly oscillatory only in its non-harmonic 
component. 

• Test Case 3: for the third test case, we work with c{y,uj) chosen ac- 
cording to dSnD, where P(cfc = 7) = P(cfc = 6) = 1/2, with 7 = 1 and 
6 = 3. The problem is thus highly oscillatory both in its non-harmonic 
and its harmonic components. 

We take the meshsize h = 0.2. The Newton algorithm is initialized with the 
solution to 

— div [{a{y, u) + c{y, co)){^ + Vwo)] = in Qn, wq is QAr-periodic, 
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and the iterations stop when 



w 



,n+l 
h 



< tol. If tol is chosen too 




large, then ( 155|) is inaccurately solved, and the variance reduction is not very 
good. For our numerical tests, we set tol = 10~^: the discrete problem ( 155|) 
is accurately solved, while only a limited number of iterations (in practice, 
around 5 iterations) are needed. 

For the numerical tests, we adopt the convention that Qn = (~^/2, 
For each Qtv, the standard Monte Carlo results have been obtained using 
2M = 100 realizations (from which we build the empirical estimator f|T3|) ). 
For the antithetic variable approach, we have also solved 2M corrector prob- 
lems, from which we build the empirical estimator f[T^ . Therefore, in all 
what follows, we compare the accuracy of the Monte Carlo approach (MC) 
and the Antithetic Variable approach (AV) at equal computational cost. 



In this test case, the energy density is positively homogeneous. We therefore 
know, from Proposition [T] and Remark [TT], that our approach yields esti- 
mations of the expectation of W^{u,C)^ i ' d^W^{uj,C) and i^d'^W'^{u,^)i 
with a smaller variance than the standard Monte Carlo approach. Our aim 
here is to quantify the efficiency gain. Note also that we have not taken into 
account, in our implementation, the fact that W^iu,^), ^ ■ d^W'^iuj,^) and 
d'^W^{(jj,^)^ are here proportional to one another. 

To begin with, we show on Figure [1] the estimation by empirical means 
(along with a 95 % confidence interval) of several homogenized quantities 
(the energy density, its derivatives with respect to each component of ^, 
. . . ) . We observe that the variance of all quantities decreases when the size 
oi Qn increases, and that confidence intervals obtained with the antithetic 
variable approach are smaller than those obtained with a standard Monte 
Carlo approach, for an equal computational cost. 

We now turn to a more quantitative analysis of the variance. Figure [2] 
shows the variances 



as a function of (note the factor 1/2 in the definition of Vyic-, consis- 
tent with f|T2l) . f|T3l) and f|T^ ). We observe that the variance of any of our 
quantities of interest (obtained either with the Monte Carlo approach or the 



3.3 Test Case 1 



Var[l^^(-,0] and V^av = Var [l^^(-, 



(56) 
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Figure 1: Test-Case 1: Homogenized quantities as a function of 
A^, for the vector ^ = (1, 1)-^ (Blue: Monte Carlo results; Red: 
Antithetic variable approach; Dashed lines: 95% confidence interval, 
equating the cost of the two approaches). From top left to bot- 
tom right: E[W^^(-,0], E[%W^^(-,0], E[d^,WU;Ol E 0], 
E[d^,^,WU;0], E[%5,W^(-,0], E[C-d^WU;0] and E[e^9|l^;y(■,0e]• 
27 



Antithetic Variable approach) decreases at the rate I/IQat] as increases 
(as expected if one could use the Central Limit Theorem). We also observe 
that the variance obtained with our approach is systematically smaller than 
the Monte Carlo variance, in the sense that Vav < Vuc- 
We next report on Table [1] the variance reduction ratio 

^AV 2Var 

which measures the gain in computational cost at equal accuracy, or the 
square of the accuracy gain at equal computational cost. Although this ratio 
somewhat varies with A^, we observe that it is of the order of 10 for all 
quantities of interest, except for d^-^^^^W^, for which it is always larger than 
4. In particular, even if A^ is not large (because we cannot afford to work on 
a large domain Qn), we still observe variance reduction. 



2N 












•966 


e ■ d^W]^ 




10 


19.41 


11.26 


13.86 


9.846 


5.966 


13.34 


19.39 


19.41 


20 


22.82 


11.89 


13.03 


9.865 


7.306 


9.096 


22.77 


22.83 


40 


18.08 


11.82 


9.816 


9.576 


5.904 


8.831 


18.03 


18.11 


60 


21.26 


12.89 


12.98 


10.57 


7.247 


10.73 


21.24 


21.28 


80 


12.36 


8.798 


9.050 


10.05 


4.316 


8.454 


12.31 


12.37 


100 


11.88 


9.856 


8.412 


11.10 


3.775 


10.24 


11.82 


11.88 


200 


13.60 


8.261 


11.52 


8.057 


4.636 


12.62 


13.54 


13.61 



Table 1: Test-Case 1: Variance reduction ratios (1571) . 



Remark 13. Similar variance reduction ratios are obtained in the case when 
the corrector problem is supplemented with homogeneous Dirichlet boundary 
conditions on the boundary on Qj\f, rather than periodic boundary conditions 
as used here following ([8]) (results not shown). 

3.4 Test Case 2 

We now consider a test-case for which the energy density is not positively 
homogeneous. From our results of Section 12. H we know that our approach 
yields variance reduction for the estimation of E [W^Ar(''0]- ^i'^ here is 
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Figure 2: Test-Case 1: Variances ( l56l) of the same quantities of interest as on 
Figure [H as a function of (Blue: Monte Carlo approach; Red: Antithetic 
Variable approach). 



29 



two-fold: we first quantify the efficiency gain, and we next verify (and tfiis 
will indeed be the case) that we also obtain a gain in efficiency for quantities 
of interest (such as the first or second derivatives of W^{u,^) with respect to 
C,) for which we do not have theoretical results in the two-dimensional case. 

We show on Figure [3] the variances fl56|) of the same quantities of interest 
as previously (obtained either with the Monte Carlo approach or the Anti- 
thetic Variable approach). As for the previous test-case, we observe that all 
variances decrease at the rate I/IQtvI as increases. In addition, we observe 
that the variance obtained with our approach is systematically smaller than 
the Monte Carlo variance, in the sense that Vav < ^mc- 

On Table [21 we report the variance reduction ratios (157]) (with the same 
convention as in Tabled]). We observe an efficiency gain of more than 10 for 
all quantities of interest, except again the cross derivative d^^^^^^^, for which 
the gain is smaller, and of the order of 4. 



2N 














e ■ d^W^ 




10 


20.38 


11.57 


14.14 


9.940 


6.206 


13.28 


19.89 


19.57 


20 


23.86 


12.34 


13.32 


9.993 


7.548 


9.265 


23.33 


23.00 


40 


18.94 


12.16 


10.16 


9.726 


6.060 


8.902 


18.50 


18.24 


60 


22.11 


13.30 


13.35 


10.73 


7.513 


10.88 


21.68 


21.41 


80 


12.89 


9.080 


9.295 


10.09 


4.420 


8.598 


12.61 


12.45 


100 


12.37 


10.17 


8.635 


11.21 


3.896 


10.24 


12.12 


11.96 



Table 2: Test-Case 2: Variance reduction ratios (157]) . 



3.5 Test Case 3 

We eventually turn to our final test-case, where both coefficients a and c do 
depend on the space variable. 

We show on Figure H] the variances (!56|) of our quantities of interest. 
Again, we observe that they all decrease at the rate I/IQatI as A^ increases, 
and that the variance obtained with our approach is systematically smaller 
than the Monte Carlo variance. 

On Table [3l we report the variance reduction ratios flSTj) (with the same 
convention as in Table [1]). Results are quantitatively similar to the ones 
obtained on Table |2) we do observe a robust variance reduction, even in 
cases for which theoretical support is still currently missing. 
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Figure 3: Test-Case 2: Variances ( 156|) as a function of (Blue: Monte Carlo 
approach; Red: Antithetic Variable approach). The quantities of interest are 
the same as on Figure [H 
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Figure 4: Test-Case 3: Variances ( 156|) as a function of (Blue: Monte Carlo 
approach; Red: Antithetic Variable approach). The quantities of interest are 
the same as on Figure [H 
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2N 














e ■ d^wt, 




10 


14.26 


12.69 


10.00 


12.38 


8.333 


10.65 


14.76 


19.37 


20 


10.82 


8.166 


7.669 


8.304 


7.730 


8.827 


11.29 


18.11 


40 


7.014 


7.077 


5.613 


10.28 


6.776 


7.310 


7.731 


14.32 


60 


10.45 


10.84 


8.666 


11.72 


8.896 


9.524 


11.82 


19.01 


80 


6.961 


5.880 


7.250 


8.800 


4.646 


8.996 


7.522 


11.10 


100 


8.543 


6.780 


7.970 


8.873 


4.669 


10.26 


8.798 


11.66 



Table 3: Test Case 3: Variance reduction ratios (!57|) . 
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